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We present an information geometric characterization of Grover's quantum search algorithm. 
First, we quantify the notion of quantum distinguishability between parametric density operators 
by means of the Wigner-Yanase quantum information metric. We then show that the quantum 
searching problem can be recast in an information geometric framework where Grover's dynamics 
is characterized by a geodesic on the manifold of the parametric density operators of pure quantum 
states constructed from the continuous approximation of the parametric quantum output state in 
Grover's algorithm. We also discuss possible deviations from Grover's algorithm within this quantum 
' information geometric setting. 
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f— - n ' Information geometry is the application of differential geometric techniques to the study of families of probabilities, 
, both classical and quantum, either parametric or nonparametric Much of the work in its quantum version 
Oh 1 has been concentrated on manifolds of density operators for both finite and infinite dimensional quantum systems. 
Classical and quantum Fisher information metrics on statistical manifolds play a key-role in information geometry, 
especially in the geometric and informational approach to classical and quantum estimation theory, respectively 

S Applications of information geometric techniques appear in statistical mechanics, in the study of the three-dimensional 
spherical model 0, [|[ , the Ising model on planar random graphs [!, H[ and the one-dimensional Potts model [!, [|| . 
Furthermore, the Fisher information plays a fundamental role in quantifying Heisenberg's uncertainty principle from a 
, statistical inference point of view [6|. Indeed, information geometric methods in combination with probable inference 
■ techniques have been recently employed to define the notion of distinguishability between probability distributions 
of curved statistical manifolds underlying the information dynamics used to provide a probabilistic description of 
physical systems when only partial knowledge (incomplete information) about them is available • This line of 

investigation is very stimulating and appealing, especially considering that it is an old dream that of viewing quantum 
mechanics as rooted in making statistical inferences based on observed experimental data [IJ Ell • Recent works where 
information geometry and inference methods are used to investigate the origin of fundamental theories as information 
geometric inferential theories appear in [T^-EBl- 

Very recently, we have observed that information theory has motivated conventional differential geometric inves- 
tigations of quantum state spaces. For instance, Riemannian geometric tools have been used to characterize the 
quantum gate complexity in quantum computing [l6-18]. In jig ], the problem of finding optimal quantum circuits 
was recast as a differential geometric problem. Specifically, a Riemannian metric on a manifold of multi-qubit unitary 
transformations is used to define a metric distance between the identity operator and the desired unitary operator 
representing the quantum computation. The metric distance turns out to be equivalent to the number of quantum 
gates needed to represent that unitary operator. It was then shown that finding optimal quantum circuits is essentially 
equivalent to finding the shortest path (geodesic) between two points in a certain curved geometry. However, the 
number of investigations concerning the applications of quantum information geometric methods to specific problems 
of quantum computing, such as quantum estimation, quantum channel entropy, quantum noise and entanglement 
and quantum algorithms is very limited. This is not unreasonable considering that the systematic study of quantum 
information geometry has only started no more than ten years ago [1] . 

In this paper, we provide a simple application of quantum information geometry to quantum computing. It is known 
that Grover's quantum search algorithm [19] can be viewed as a geodesic path on the manifold of Hilbert-space rays 
where the notion of quantum distinguishability is quantified via the Fubini-Study metric, a gauge invariant metric 
on the projective Hilbert space [2(], Hl| . By observing that a parametric quantum wavefunction induces in a natural 
manner a parametric density operator and by considering its square root, we use the concept of Wigner-Yanase 
quantum information metric. Such metric is one among many versions of a so-called quantum Fisher information 
metric, a metric on manifolds of density operators for both finite and infinite dimensional quantum systems. We 
show in an explicit manner that the Wigner-Yanase metric and the Fubini-Study metric differ by a factor of four 
when considering pure state models. Finally, interpreting the Fubini-Study metric as a quantum version of Fisher 
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metric, we provide an information geometric characterization of Grover's algorithm as a geodesic (shortest length 
curve) in the parameter space characterizing the pure state model, the manifold of the parametric density operators 
of pure quantum states. Finally, we discuss few possible deviations (non-constant Fisher information function and/or 
non- actuality) from Grover's algorithm within this quantum information geometric framework. 

The layout of the article is as follows. In Section II, we introduce the Wigner-Yanase quantum information metric 
from both formal and heuristic viewpoints. In Section III, we briefly discuss Grover's search algorithm in quantum 
computing. In Section IV, we show that Grover's dynamics corresponds to the shortest path from an information 
geometric point of view. In Section V, we consider deviations from Grover's probability path where the requirements 
of constant Fisher information function and/or actuality are not simultaneously satisfied. Our concluding remarks 
appear in Section VI. 



II. ON INFORMATION GEOMETRY AND STATISTICAL DISTINGUISHABILITY 

In the classical setting, it is known that except for an overall multiplicative constant the classical Fisher information 
metric is unique p2l |23| : it is the only monotone Riemannian metric with the property of having its line element 
reduced under Markov morphisms (stochastic maps). Said otherwise, there is essentially one classical statistical 
distance quantifying the classical distinguishability between two probability distributions. In the quantum setting, 
Riemannian metrics are considered on the space of density matrices. The requirement that the distance between 
density matrices expresses quantum statistical distinguishability implies that this distance must decrease under coarse- 
graining (stochastic maps). Unlike the classical case, it turns out that there are infinitely many monotone Riemannian 
metrics satisfying this requirement |24j-[26|. 

In this section we introduce the Wigner-Yanase quantum information metric from both a formal and heuristic point 
of view. 



A. The Wigner-Yanase Quantum Information Metric: the Formal Approach 

Assume that M. n denotes the space of (n x n)-complex matrices, T> n is the manifold of strictly positive elements 
of A4 n and T>\ C T>„ is the submanifold of density matrices, 

V\ d = {D £ T> n ■ D = Z^and Tr (D) = 1} . (1) 
The tangent space to T> n at p is given by, 

T p V n d = {D £ M n ■ D = D f and Tr (£>) = 0} . (2) 
Given D\ and D 2 in T p P^, the Wigner-Yanase quantum information metric is defined as [27j . 

(D lt D 2 ) p: /wy d ^ Tr [Dj • c WY {L„, R p ) (D 2 )} , (3) 
where L p and R p are the left and right multiplication operators, respectively, with 

L p (D) d = pD and, R p (D) = Dp. (4) 

The quantity cwy (L p , R p ) in (|3]) denotes the so-called Chentsov-Morozova function for the Wigner-Yanase informa- 
tion metric, 

cwy (x, y) = f 1 . s = -2 with x, y > 0. (5) 

y/wv (J J W% + y/y) 

The function / : (0, + oo) — > M is an operator monotone function, that is Vrt € N and V Mi, M 2 £ M. n with 
< Mi < M 2 , it must be < / (Mi) < /(M 2 ). To evaluate (D 1 , D 2 ) /wy in ©, observe that T p V x n can be 

decomposed in two orthogonal subspaces (T p V\^ c and (T p I?^) , 

rX = (rX) c ©(T p pi)°, (6) 
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where, 



{T p Vl) c = {De T p Vl : [D, p] = 0} 



(7) 



and (TpX>*)° is the orthogonal complement of (T P X>^) C . Therefore, an arbitrary element I?' of (TpT>]^ a can be written 
as i [p, D] with D = D^. Thus, 



t P , /wv = (* ^ > * l>- °]) Pl / WY = Tr \P> °] 4 ( x p + R ?) ( l I/ 3 ' ^) 
where we have made use of the following two relations, 

[p, D] = (L p - R p ) (D) = (4 + Rj) (if - 4) (D) , 



= -4Tr 



D 



and, 



We notice that ID', D' 



P, fv 



L} (D) =p*D, m {D)=Dp*. 
in (|5J) can be written as, 

Jwy(p) = (D', D') ptfwY =8I(p, D), 



where I (p, D) defined as, 



P 5 , D 



(8) 

(9) 
(10) 
(11) 
(12) 



is the skew information introduced by Wigner and Yanase when studying quantum measurement theory from an 
information-theoretic viewpoint [28]. 



B. The Wigner- Yanase Quantum Information Metric: the Heuristic Approach 

The classical Fisher information lp (pg) of a parametric probability density pg (x) is a key-quantity in the statistical 
estimation theory and it is defined as, 

If (pe) = J dxp (x) \ ^ J . (13) 

Observe that If {pg) can be rewritten in terms of a fundamental quantity in quantum physics, the probability ampli- 
tude (wavefunction) y/pg (x), 

I,<J»)=*Jd*(y^. (14) 

When going from classical to quantum theory, a quantum version of the classical Fisher information can be heuristically 
defined as follows: the integral in (|14[) is replaced by the trace and the probability densities pg by density operators 
p g . Therefore lp (pg), a quantum extension of lp (pg) to a family of parametric quantum operators p e , reads 

MPe)=4Tr(^) 2 . (15) 

Quantum generalizations of the very same classical Fisher information are not unique. Indeed, two classically identical 
expressions (see (fT3"|) and (O, for instance) generally differ when they are extended to a quantum setting. This 
difference is a manifestation of the non-commutative nature of quantum mechanics and is reminiscent of the idea of 
quantum discord [29]. 

If pg in (|15[) satisfies the von Neumann equation, 



(16) 
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with p Q = p, 9 S R is a temporal parameter and T is the generator of temporal shift. Then, 



Pe = e-^pe^ and, ^ = ie~ iST [pi, T \ . (17) 

Substituting the second equation in (JTTJ) into (Tl~5l) . (p e ) becomes independent of 9 and reads 

I F (p e ) = 8I(p,T), (18) 

where / (p, T) is the Wigner-Yanase skew information defined in (|12[) . Thus, the skew information is a particular kind 
of quantum Fisher information Ip (pg) known as the Wigner-Yanase quantum Fisher information 

I WY (p e )=I F (p e ) = 8I(p,T). (19) 

When a density operator is parametrized by n-real parameters 9 = (#,..., 0") and the quantum state p e satisfies 
the von Neumann equation, the set of quantum states p g forms the quantum evolution submanifold of X>*. The set 
of parameters 9 can be regarded as a local coordinate system on such submanifold endowed with a metric structure 
defined by the information metric g^j^ W> 



4 W Y) (9) = 4Tr [(d^) (d iy /to)] , 5, = A (20) 



Once the monotone Riemannian metric g\J^ (^) i s explicitly obtained, the other geometric quantities (Christoffel 
connection coefficients, scalar and sectional curvatures, Riemannian curvature tensor, etc.) can be calculated from it 
in a purely mathematical way, at least in principle. For instance, the geodesic equation reads [30j . 

n k del dg >* 



where D denotes the covariant derivative, 9 = and, 

It (9) = \g kl (d i9lj + d m - d l9ij ) , (22) 
arc the Christoffel connection coefficients. 



III. ON GROVER'S SEARCH ALGORITHM 



In this section we briefly discuss what is considered a masterpiece of quantum computational software [3l[, that is 
Grover's search algorithm in quantum computing. 



A. Quantum Searching 



The search problem may be stated as follows: we wish to retrieve a certain item satisfying a given condition assuming 
that it belongs to an unsorted database (oracle) containing N = 2™ items. One step is needed to specify whether or 
not the examined item is the one fulfilling the given condition. We assume that the selection of the item is aided by no 
sorting on the database. In such a situation, the most efficient classical algorithm to tackle the search problem requires 
the examination of the items in the database one by one. Therefore, by means of a classical computer, the oracle 
must be queried on average ^ times ( O (N) classical steps). However, by using the same amount of hardware as in 
the classical case but by having the input and output in superpositions of states, Grover has developed a quantum 

mechanical algorithm capable of solving this search problem in about j \fN steps ( O (v^V^j quantum mechanical 

steps) fl9| . Although this is not dramatic as the exponential quantum advantage achieved by Shor's algorithm for 
factoring, the extremely wide applicability of searching problems makes Grover's algorithm interesting and important. 
In particular, Grover's algorithm gives a quadratic speed-up in the solution of ./VP-complete problems, which account 
for many of the important hard problems in computer science. Indeed, it was shown in [32j that relative to an oracle 
chosen uniformly at random, with probability 1, the class NP cannot be solved on a quantum Turing machine in time 
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O yy/Nj. Drawing on this result, Grover pointed out in [33| that his algorithm is optimal (i.e., the fastest), up to a 

multiplicative constant factor, among all possible quantum algorithms. A detailed proof of Grover's statement about 
the optimality of his quantum searching algorithm appears in [34j . Finally, in [351 ] , it is shown that for any number 
of oracle lookups, Grover's algorithm is exactly (and not just asymptotically) optimal. 

What makes a quantum search algorithm more efficient than another? An algorithm is an abstract mathematical 
concept, whereas it is useful to consider how efficiently we can run an algorithm on a computer. Computer scientists 
associate a cost with each step of the algorithm and with the amount of memory required, embodying the idea that 
physical computers have a finite size (memory) and work at a finite rate of elementary calculation steps per unit time. 
This gives us a way to determine if one algorithm is intrinsically faster than another. In general, it is very difficult to 
formally prove that a certain algorithm is the "best" algorithm for a given computational task. For instance, Shor's 
factoring algorithm [36j j is only the "best known" algorithm for factoring; there is no proof that something faster 
cannot be found in the future. 



B. Grover's Algorithm 



In what follows, we outline the construction of Grover's algorithm when considering the n-qubit case (N = 2™ 
states). The step-0 (the initialization) of Grover's algorithm begins by using the Hadamard transform to construct 
a uniform amplitude initial state which is an equal superposition of all the orthonormal computational basis states 
{|s)} in the iV-dimensional Hilbert space, 



N-l 



(23) 



Observe that the state (|23p may be rewritten as 

N-l 

N 



N-l 



N-l 

N 



F L TEl s ) + ^l«)= cos (f) |r>+ain(f)|a>, 



s— s^a 

where the angle a characterizes the overlap between \q) and the searched state |a), 



(24) 



(25) 



while the state \r) is defined as, 



V) = 



N 



(26) 



Given the initial input, a single-step of Grover's algorithm is characterized by a rotation by a in the two-dimensional 
space spanned by \r) and \a). The effect of the algorithm after m-steps leads to the state \i/j g (m)), 













^ 1 


COS 


(^m - 




\r) + sin 


- 












-2) a _ 



(27) 



In the limit of N ^> 1, the number of steps m for which \tp G (m)) coincides with the searched state \a) (with success 
probability equal to 1) is approximately given by, 



- JV>>1 n HTt 
m ~ — v I V . 
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(28) 



Equation (|2"5)l can be obtained by imposing that (m + I) a = 5 and by observing that when N ^> 1, Eq. (|2~5"j) implies 
that a ~^=. 

Although Grover's algorithm evolves with discrete m, in the limit of N ^> 1 the output state (f2"7T) can be approx- 
imated by a quantum wave-state \ip (8)} depending on a continuous parameter 9. Indeed, considering the following 
formal substitutions, 



1 



N-l 



m+ - a -t 6, \a) -> |0) , \r) 



VN 



k=i 



(29) 
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the state (9)) reads, 

N-l 

\i>m = E vWwi*). (30) 

fc=0 

where, 

{k\k') = 5 kk , , po (9) = sin 2 and, Pj (0) d ^ f |^ with j ^ 0. (31) 

The ./V-dimensional probability distribution vector p = (p (9) , p\ (9) pn-i (9)) with p k {9) defined in (|3Tj) can 
be regarded as a path characterizing Grover's algorithm on a suitable probability space. In what follows, we show 
that such a path is indeed a geodesic path for which the quantum Fisher information action functional achieves an 
extremal value. 



IV. INFORMATION GEOMETRY AND GROVER'S ALGORITHM 



In this section, we show that Grover's dynamics leading to (|30[) corresponds to the shortest path from an information 
{cometric point of view. 



A. The Parametric Pure State Model 



Consider the parametric density operator p g constructed using the quantum state \tp (9)} = \%p g ) in (I3U1) . 

N—i 

Pe = m (1>e\ = E V)Pi (°t 2 e-A^-'t'M \k) (j\ , 

k, j=0 



(32) 



where the pure state \ip g ) is normalized to one. Since p g is a pure state, it satisfies the following relations, 

Pe = Pe and > Pe = y/Pe- (33) 



We assume that the notion of distinguishability on the manifold of density operators defined in (|32[) is provided by 
the Wigner-Yanase quantum Fisher metric defined as, 



l WY 



(p g ) d = f 4Tr [(3 iV 7£) (d Jy /pi)] = 4Tr [(d lPg ) (d jPe )} 



(34) 



where in the multi-parametric case 9 = (d 1 ,..., 6 l J and di = f -J— with i = 1,..., /. The second equality in AMI) holds 
because p g is a pure state in the case being considered. After some algebraic manipulations and observing that the 



normalization condition for \ip g ) implies that, 

(djljlgllpg) = - (ipgldjtpg) , 

the Wigner-Yanase quantum Fisher metric becomes, 

(Pe) = 4 [{dii>e\dji>e) + (d^el^e) (fy^el^e) 



y-(quantum) 
1 WY 



Thus, the infinitesimal line element on the manifold of density operators in (1321) reads, 



T(quantum) 
WY 



( Pg ) d9 l d9 J = 4 [{d^d^g) + (MeWe) {d^ g \^ g )\ d9 t d9, 



(35) 



(36) 



(37) 



Defining 7, ■ and aij as, 



lii =* Re[(dii>e\djiPe)} and a l3 d = Im [{diip 6 \dj^ e )} , 



(38) 
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respectively, it follows that 



Because of the asymmetry of Cj,-, we have 



that is, 



lij = Jji and, a l} = -o yi . 
OijdOidOj = —OjidOidQj = —<7jid9jd0i — —OijdOidOj, 
dijdOidOj = 0. 



Therefore, it turns out that ds^y m (03 becomes 



d<? 

ui >WY 



(quantum) 
WY 



( Pe ) dOidOj = 4 [Re (Wold^e) + (ftVshM (drfeWe)] d9 l d9 



We remark that ds^y m (1421) is exactly four times the Fubini-Study infinitesimal line element is given by 37 [ 



ds 



FS 



g p (9)de i d9 3 = iidv.ii 2 - kvi#)| 2 = i - K^'iv-)! 2 , 



where \ip) = \ip (9)) = \ip d ), \4>') = \4> (9 + d9)) = jV'e+de) an d the Fubini-Study metric g^ S ^ (9) reads 

5 jf s) (9) ee 9ij (9) d ^ [Re (ftVflia^*) + (SflMlM (d^eWe)] ■ 



(39) 
(40) 
(41) 

(42) 
(43) 
(44) 



The Fubini-Study metric is a gauge invariant metric on the manifold of Hilbert-space rays (projective Hilbert space). 
Our analysis explicitly recognizes that the Fubini-Study metric is a quantum version of the Fisher metric |38| . 



B. Grover's Information Geometric Dynamics 

Substituting (|4"3"f into (l4"2l and using (|30[) together with the normalization condition on the probabilities in (|31l) . 
after some straightforward algebra it turns out that the infinitesimal Wigner-Yanase line element reads (for further 
details, see Appendix A), 



tts WY 




N-l 



(N-l \ 
k=0 J 



d9 z 



where. 



. dof dp k (9) 

Pk = and 



def #fc (6) 

d9 ' 



(45) 



(46) 



Recall that Grover's algorithm consists only of a sequence of unitary transformations on a pure state. In the continuous 
approximation, such sequence of quantum unitary operations leads to the output state \4> {9)} in ()30p . Such quantum 
mechanical wave- vector is such that its Fisher information function J- (9) [39( is independent of the parameter 9 (it 
is a constant quantity), 



N-l / ~i \2 N-l . 2 N-l /a , N 2 

k=0 V ' k=0 F k=0 



89 



4. 



(47) 



We point out that the general expression of the Fisher information function T (9) in (|47[) is invariant under unitary 
transformations. This can be explained as follows. Consider a normalized pure state \ip (9)} given by, 



N-l 



foH0)> = E VM9)e' 



<t> k {8) 



\k) 



k=0 



Eq. (gSJ) implies that, 



IV>0 ( m )l = VPm (9), 



(48) 



(49) 
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dcf 

where i/jq (m) = (mlV 7 (0)) with (m\k) — 5 m k and the Fisher information function reads, 

m . t 'f(?l*^)'. K f&. (50) 

m=0 v 7 m=0 ^ 

Under unitary transformations J7 mapping \tp (9)) to (0)) = U \ip (0)), the transformed Fisher information function 
F' (9) becomes, 

m— \ / m— 

where i/'e ( u m) == (u m \ip' (9)) with (u m |ufe) = 5 m k and = f t/ \k). Thus, from (|50|) and (jSTj) it turns out that the 
Fisher information function remains unchanged under unitary operations. In the rest of the manuscript, we refer to 
Eq. (|47j) as the "parametric-independence constraint" on the Fisher information function. 

As a side remark, we emphasize that within a statistical inference viewpoint extended to the quantum framework, 
the statistical notion of Fisher information function resembles a generalized mechanical notion of kinetic energy with 
respect to a given statistical parameter (0, in our case) regarded as temporal or spatial shift. Specifically, it can be 
shown that the following relation holds [40], 

K{e) = \T{e) + Y d J 2 e {k)\^ e {k)\\ (52) 

k=0 

with T (6) defined as in (I5TJ1) and ip e (fc) = (k\ipg) and p k (0) = \ipg (k)\ 2 . The quantity JC (9) denotes the kinetic 
energy of the wavefunction with respect to the parameter 9 while Jg (k) is a statistical analogue of the normalized 
quantum mechanical current density with respect to 9. They are defined as [40], 



jc (9) y: 



def 

k=0 

and, 



^ 1 ^ e (fc)| 2 



00 



(53) 



respectively. The symbol "*" denotes complex conjugation. For ^ g (k) = yjp k (9)e^ k< ^ e \ we obtain that 

jv-i . 2 

fC (9) = (^ e \4, e ) ,T(9)=Y.- and > Je (*0 = k (V) > (55) 

k=0 Pk 

with %jj e = f ■ Therefore, Eq. (f5"2"j) reads 



do 

N-l .9 N-l 



4 ^ Pk 

k=o 1 fc=0 

Since in Grover's case 7 s (9) = 4 and Jg (k) = for any k — 0,..., N — 1, condition ([56]) becomes 

K{0) = (i> e \^ e ) = l. (57) 



We conclude that from a statistical quantum inference viewpoint, Grover's algorithm is characterized by a constant 
statistical kinetic energy where no "statistical dissipation" occurs. 

The geodesic path related to Grover's algorithm is found by minimizing the action S [p k (9)] defined as, 



S [Pk (9)} = / a/^wy = / £ (Pk (0) , Pk (9)) d9, (58) 



9 



where the Lagrangian-like quantity C (p k (9) , Pk (9)) reads, 



C{ Pk (0),p k (9)) 



£ 

k=0 



vl (0) 



Pk (0) 



given the normalization constraint on the parametric probabilities pk (9) , 



N-l 



(59) 



(60) 



k=0 



For the sake of simplicity, consider the change of variable pk {9) — > qi (9) (4l|. Then, using the Lagrange multipliers 
method, the new action S' [qk (9)} to minimize becomes , 



S'{q k {9)\ = / C , (q k (9),q k (9))d9 = 



N 



L k=l 



N 



-A £g2(tf)-l 



d<9, 



(61) 



where A is the Lagrange multiplier and the new Lagrangian-like quantity is given by, 



C (q k (9) , q k (9)) 



N 



AT 



-A £^(0)-l 



vfc=l 



The path minimizing the action S' [q k (9)] satisfies the " actuality constraint" , 

SS' [q k (9)] 



5q k (9) 



= 0, 



leading to the Euler-Lagrange (EL) equation, 



d_ 

d0 



d£'(q k (6),q k (d))\ dC'(q k (9),q k (9)) 



dq k 



dq k 



From (|fj2"j) . it follows that 



dC 
d9 \dq k 



4q k 



qkqk 



1 



N 

4E€ 

k=l 



and, 



dC 

dq k 



N 

fc=i 



-2\q k . 



N 

fc=i 



= 0. 



■4g fc , 



(62) 



(63) 



(64) 



(65) 



Thus, substituting (|66[) and (|65j) into (|64l) . the EL equation becomes 

d 2 q k {9) l{q k {9),q k {9))dq k {9) A 



~£fe(0),<ft(0))<ft(0) = O, 



C{q k {9),q k {9)) d9 

_ dC{g k (0), q k (9)) „ rUU rfA n _ l0 ^, ,.;,.,„, ; „ , r ,,n „.;,!, ,,'-! 

that £ (^fc (0) , q k (9)) is constant and equals two while C (q k (9) , q k (9)) equals zero. It finally follows that 

cos 9 



(66) 



(67) 



where C (q k (9) , q k (9)) = aL -^°J g 3hm> w j t h £ (q k (9) , q k (9)) given in ((5SJ with g 2 (0) = p fc (9). Using g7|), it follows 



go (9) = sin and, q k (9) — 



(68) 



where k = 1,..., N — 1 fulfil the Eq. (|6"T|) with A equal to one in order to satisfy the normalization condition (f6"0|). We 
then conclude that the iV-dimensional probability vector p = (po (6*) , pi (0) pn-i (9)) with p/j (9) defined in (|3"Tj) 
is a geodesic path for which the quantum Fisher information action functional S' [q k (9)] achieves an extremal value. 
We also stress that such a path satisfies both the parametric-independence and actuality constraints in Eqs. (1471) and 
(f6"3")l . respectively. 
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V. DEVIATIONS FROM GROVER'S ALGORITHM 



As pointed out earlier, Grover's search algorithm relies on two key-assumptions: 1) Grover's iteration must be 
characterized only by means of unitary transformations acting on a uniform amplitude input state. This leads to 
the parametric- independence constraint on the Fisher information function (see Eq. (|47p ); 2) the expressions for the 
components of the probability distribution vectors associated to Grover's algorithm are such that the quantum Fisher 
information action functional is minimized when evaluated on such paths (see Eq. (|63l0 . In summary, Grover's paths 
are actual unitary probability paths of constant Fisher information. When one of these two requirements (parametric- 
invariance and actuality) do not hold anymore, we are led to "deviations" from Grover's algorithm which should be 
somehow less efficient than Grover's one. 

In what follows, we discuss unitary deviations from Grover's algorithm and attempt to quantify, whenever possible, 
different aspects from an information geometric point of view on manifolds where the statistical parameter 9 can be 
regarded as a local coordinate parameterizing the density operator p g . 



A. Grover's Model: Actuality and Parametric-Independence 



Recall that the probability distribution vectors associated with Grover's algorithm is given by, 

9 „ cos 2 9 cos 2 i 



p G = sin ti, 



Eq. (loT)|) leads to the constancy of the Fisher information function F (9), 

= 4. 



k=0 



Pk V dO 



The geodesic equation in (f2~Tj) for Grover's statistical model reads, 

2 , i *rn? = Qf 



IT (9) d9 



(69) 



(70) 



(71) 



with 9 



del d6(r) 
dr 



Since J- (9) is constant, the geodesic equation reads 9 = 0. Assuming that the boundary conditions 



are given by 9 (n) = 9i and 9 (t/) = 9f, integration of (fTTj) leads to the geodesic path 9q (t), 



0g (r) 



iTf - tip 



(72) 



Tf ~ T t 

The geodesic path 9q (t) can be regarded as a continuous succession of intermediate states connecting 9i to 9f. The 
geodesic motion occurs with speed vq (t), 



"Grovor ( T ) 



del' 



d9 Gl 



dr 



satisfying the normalization relation, 



9lm {S) 



d9 l (t) d9 m (r) 
dr dr 



= 1, 



(73) 



(74) 



with I = m = 1 and gn (9) = 4. Equation (|74]) implies that the temporal duration At needed to navigate the geodesic 
path connecting 9i to 9f is given by, 



def 



Ar = T f -T i = 2(e f -e i ). (75) 
We point out that At can also be viewed as the length Lq (9i, 9f) of the geodesic path connecting 9i to 9p 



del 



L G (6i,e f )= I gim (9)d9 l d9 m 
Indeed, in the case being considered Lq (9i, 9f) reads 



dr 



9lm {0) 



d9 l (t) d9 m (t) 
dr dr 



L 



dof 



G<»z, Of) 



Vld9 = 2(9 f - 9i) = At. 



(76) 



(77) 
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B. Deviation from Grover's Model: Neither Actuality Nor Parametric-Independence 



In this example, we consider a properly normalized probability distribution vector given by, 



PModcl-II 



1 - 9 2 1-9 2 



N-V N - 1 



(78) 



Notice that when 9 <C 1, pbrover ~ PModel-ii- The Fisher information function T (9) for the modified Grover's statistical 
model is not constant (no parametric-independence) and reads, 



i fd Pk (ey 2 



k=0 



Pk V dQ 



(79) 



We also remark that PModel-n with T (9) in (1791 does not satisfy the EL equation in ([571) . According to the principle 
of least action, this means that p*Modei-ii can be regarded as a possible but not an actual (not extremal) probability 
path. Furthermore, the geodesic equation ([2~Tj) for the modified Grover's statistical model becomes, 



= 0, 



(80) 



dr 



where 8 

(I5U1) becomes a first order differential equation, 



Eq. (|80l) is a nonlinear second order ordinary differential equation. Setting 9 = —^p- = y (0), Eq. 



1)1) 



tV 2 = 0, 



(81) 



where y' — dv j^ ■ Assuming y ^ 0, integration of (|5T1) implies that 



y{9) = C{\-9 2 )\ 



(82) 



where C € K is an integration constant. Since y (9) — de } T ' , it finally turns out that the general solution of Eq. ([80 
is given by, 



'(r)=sin [C (r + Ci)] 



(83) 



where Co and C\ are real integration constants whose explicit expressions can be obtained once the boundary condi- 
tions are fixed. Assuming 9 (t^) = 9i and 9 (r/) = 9f, the expression for the geodesic path 6*Modci-ii (t) reads, 



#Modd-ii (t ) = sin 



arcsm t) t — arcsm ( 



r f arcsm t/.; — r\ arcsm a f 
arcsin Or — arcsin 



(84) 



The geodesic path 9 (r) can be regarded as a continuous succession of intermediate states connecting 9i to 9 f . The 
geodesic motion occurs with speed v Modei-ii (t) , 



fModel-II {t ) 



dof rf^Modcl-II (t ) 



satisfying the normalization relation, 



9lm 



d9 l (t) d9 m (r) 
dr dr 



= 1, 



(85) 



(86) 



with I = m = 1 and gu (8) = 4 (l — 6* 2 ) . Equation ([5o]l implies that the temporal duration At' needed to navigate 
the geodesic path connecting 8i to 9f is given by, 



At' = f t'j — t ■ = 2 (arcsin 6*/ — arcsin ^j) . 



(87) 
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We point out that At' can also be viewed as the length L Of) of the geodesic path connecting 9i to 9f, 

d9 l (r) d9 m ( T )~ 



L 



Modol-II Vf 



del' 



9lm (6) de l de r 



dr 



9lm ifi) 



dr 



dr 



Indeed, in the case being considered I/Modei-ii {9i, Of) reads 



-^Model-II {9i, Of) 



def 



1 



—d0 = 2 (arcsin(9/ — arcsin#i) = At 



(89) 



Comparing Grover's model to the deviation from Grover's model and assuming 6f = 9i + e with 9i = and < e « 1, 
we obtain 



ATModcl-II - ATGrover + ° i^) ^ °- 



(90) 



We may conclude that navigating along Grover's geodesic path is more efficient than navigating along the geodesic 
path characterizing the modified statistical model. Indeed, assuming that any navigation that is being compared 
occurs with the same unit speed, the efficiency of the navigation is simply quantified in terms of the length of the 
geodesic path connecting the initial and final states. Furthermore, since the length of a geodesic path can be regarded 
as the navigation duration, the most efficient unit-speed navigation is characterized by the shortest navigation duration 
(shortest geodesic path) for a given set of boundary conditions. 



C. Deviation from Grover's Model: Parametric-Independence without Actuality 

It is possible to consider normalized probability distribution vectors PModel-m = (Po (0) , p\ (0) Pn-i (0)) that 
satisfy the parametric-independence constraint but do not minimize the action functional S' [qu (9)]. For instance, 
assuming that pk (9) — pv (9) for k and k' 6 {1,..., N — 1}, the components of PModci-in should satisfy the following 
system of two equations, 

Po (0) + (N- l) P - k (9) = 1 and, ^ + (N - 1) ^ = 4, (91) 
with k = 1,..., N — 1. From the first equation in (|9"Tj) , we get 

p- k (0)= 1 -j^. (92) 

Substituting into the second equation in (liJl"l) . we obtain 

J j2_4p (l-Po) = 0. (93) 

Integrating and using (|9"2")l , a solution of the system (EJll) is given by, 

_/l + sm26> l-sin26> l-sin26»\ 
PModei-m - g , 2{N _ iy -: 2{N _ 1) ) ■ ( 94 ) 

In this case, the Fisher information function J- (9) = 4, however the components of PModci-iii do not satisfy the EL 
equation in (|67|) . Although the integration of the geodesic equation in (l2lT) for this model would lead to geodesic 
paths that are straight line trajectories, the components of PModci-in represent only possible and not actual trajectories 
of the system (according to the principle of least action). Thus, the probability path pModci-in lacks the actuality 
requirement that distinguishes Grover's probability path. Indeed, a part from an irrelevant constant factor, it turns out 
that the probability distribution vector related to Grover's algorithm is the only one satisfying both the parametric- 
independence and minimization condition (actuality constraint) simultaneously. 
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D. Deviation from Grover's Model: Actuality without Parametric-Independence 



For the sake of completeness, we mention one more possible scenario. There might exist a probability distribution 
vector PModci-iv = (po (9) > Pi (9) ,—, Pn-i (9)) with p k (9) = p k > (9) for k and k' € {1,..., N — 1} that satisfies the 
minimization without fulfilling the parametric-independence requirement. The components of such a path should 
satisfy the following system of nonlinear ordinary differential equations, 

?o (9) - §§go (0) + ^q (9) = 0, 

(95) 

where q 2 (0) = po (9) and q? (9) — p^ (9) for k = 1,..., AT — 1 with the normalization condition 

ql{9) + {N-l) q l{9) = l. (96) 

The quantity £ (9) in (l95j) is defined as, 



dcf , 



£(<?) = 2 ?o (^) + (AT — 1) q (9) ^constant. (97) 



Using the normalization condition (|96|) . the expressions for £ (6») and t (9) d = ^ read 



C(9) = 2q , 



' (jV-l)(l-gg) + gg 
(iV-l)(l- 9o 2 ) 



and 



- 1) (1 ~ go 2 ) + go 2 | 2g gg (JV-l)(l- g g) 

£ {6) ~ ""I (N-l)(l-qi) + (N-l)(l- q tf V (^~D(l-go 2 ) + go 2 ' (99) 

respectively. It is fairly clear that finding an analytical expression for actual TV-dimensional PModci-iv m the absence 
of parametric- independence may turn out to be quite nontrivial (for further details, see Appendix B). It is interesting 
to notice that the same computational challenge occurs when attempting to find analytical solutions of a non-unitary 
quantum walk 42]. However, focusing our attention on probability paths p with only two types of components, po 
and p k with k = 1,..., N — 1 and A^ = 2, we are able to show that the only possible probability distribution vector 
satisfying the system (1951) is Grover's path where C (9) is constant (and equals two). In such scenario, Eqs. (|98[) and 
(j9"9]l become 

C (9) = *t and, t (9) = ^ (l = ^ + 3 ^ . (100) 

v/i^f (i-g 2 ) 3 

Using (|100p together with the first equation of (|9"5)) , we get 



go- ^l-g 2 = 0, (101) 

that is, 

<?o (0) = sin {9 + c ) and q x {9) = cos (0 + c ) (102) 

where Co is the real integration constant. We then conclude that in this special toy-case (one-qubit quantum Hilbert 
space, %\ with dimc?^ = 2 = N), actual paths are necessarily parametric-independent unitary paths as well. 
However, the realistic scenario implies a continuous approximation of Grover's quantum searching problem where 
N 1. In such a case, actual parametric-dependent probability paths with the two-component structure may still 
exist (see Appendix B). 

As a final remark, we point out that our information geometric analysis may be of use to investigate the possibility 
of the existence of actual parametric-dependent probability paths with a multi-component structure (more that two 
types of components). Such scenarios occur in a multi-item quantum search where more than a single item is being 
searched. 
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VI. CONCLUDING REMARKS AND OPEN ISSUES 



In this article, we have quantified the notion of quantum distinguishability between parametric density operators 
by means of the Wigner-Yanase quantum information metric. We then presented an information geometric charac- 
terization of Grover's quantum search algorithm as a geodesic in the parameter space characterizing the pure state 
model, the manifold of the parametric density operators of pure quantum states constructed from the continuous 
approximation of the parametric quantum output state in Grover's algorithm. Finally, we discussed few possible de- 
viations from Grover's algorithm within this quantum information geometric framework. For instance, in the second 
example in Section V, we have shown that navigating along Grover's geodesic path is more efficient than navigating 
along the geodesic path characterizing the modified statistical model. This statement holds when it is assumed that 
any navigation that is being compared occurs with the same unit speed and, therefore, the efficiency of the navigation 
is simply quantified in terms of the length of the geodesic path connecting the initial and final states. Furthermore, 
since the length of a geodesic path can be regarded as the navigation duration, we have also pointed out that the 
most efficient unit-speed navigation is characterized by the shortest navigation duration (shortest geodesic path) for 
a given set of boundary conditions. 

The quantum information geometric techniques employed in our work were fairly simple. Unfortunately, we stress 
that there exist some technical difficulties in handling quantum information geometric techniques in general cases. 
For instance, there is no general formula for geodesic paths of manifolds of density operators endowed with an 
arbitrary given monotone Riemannian metric |43| . Explicit expressions are known only for two special metrics: the 
Bures metric [44[ and the Wigner-Yanase metric (27l |. In [13], geodesic paths on manifolds of density operators are 
obtained as projections of large circles on a large sphere within the purifying Hilbert-Schmidt space. In [27| . geodesic 
paths are obtained exploiting the classical pull-back approach to Fisher information metric via sphere geometry. If 
diagonalization procedures are used, it is possible to obtain an explicit expression for the Bures distance of two 
arbitrary density matrices in the two-dimensional case (45j . In the general n-dimensional case this is unfortunately 
not possible. Thus, one of the greatest mathematical problems to face in order to quantify the metric properties 
of arbitrary mixed states is that of the explicit computation of eigenvalues of arbitrary Hermitian density matrices. 
The scenario becomes even more difficult if we recall that in order to find the geodesies one needs to integrate the 
geodesic equations expressed in terms of the metric and the Christoffel connection coefficients (which are expressed in 
terms of the first order changes in the metric). Furthermore, to integrate the Jacobi-Levi-Civita equation of geodesic 
spread, one also needs the explicit expression of the Riemannian curvature tensor (which encodes information about 
the second order changes in the metric tensor on the manifold of density operators). It may be that the first and 
second order perturbation theory of linear operators may turn out to be useful in estimating the changes in the metric 
tensor. However, it appears that thus far this line of investigation has not provided useful results |45j |. There have 
been other attempts to compute in finite dimensions explicit formulae for the Bures metric where no diagonalization 
procedure is employed. These approaches only make use of the theory of matrix equations, determinants and traces 
46] . Needless to say that further research is needed to fully handle the quantum information geometric formalism. 

As pointed out in [l|, quantum information geometry is only in its infancy, and much more research awaits to be 
performed. However, from the preliminary results obtained in this article, we have reason to think that quantum 
information geometric techniques may turn out to be especially useful for describing and understanding the efficiency 
of both unitary and non-unitary quantum walks used to implement quantum search algorithms in quantum computing 
[47j . In particular, our analysis opens up new lines of investigation that may deserve some attention. For instance, 
it would be worth deepening our analysis and attempting to understand the connection between computational 
complexity classes of quantum algorithms [48| and the complexity of the quantum geodesic paths associated with 
them. These investigations are left for future works. 
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Appendix A: Derivation of the Wigner-Yanase line element 

Here we compute in an explicit manner the infinitesimal Wigner-Yanase line element ds^yyi 

ds^y = 4ds| S > ( A1 ) 
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where the Fubini-Study infinitesimal line element ds| g is given by, 

d4 s = ll#il 2 -l^i#>i 2 = i-|(v , l^>| 



(A2) 



where \ip) and Ji/'') are neighboring normalized pure states expanded in an orthonormal basis {|fc)} with k € {!,•••, N}, 



N 



N 



J2VPk(0)e l4, ^\k) and, \^') = VPk + d Pk e 



i{<f> k +d<l> k ) 



\k) 



fc=l 



k=l 



respectively. Observe that, 

N 



N 



W) = E VPk+d Pk e l ^ +d ^ \k) = 



k=l 
N 

E 

fc=i 



fe=i 



_ ldp k I dp 2 
2 Pk 8 p£ 



|fc). 



I*) 



Eq. (|A4|) implies that, 
and thus \ip) becomes, 



N r 



fc=l 



__ 1 dp k 1 dp 2 

IPk 1 + o Q~ 

2 Pfc 8? 



±) e-^ (l - id<j> k - \d<j>l 



AT 



Pfc U + 



1 dp fe 1 dp 2 

2 8^1" 



1 — ic 



>k ~ j ^ fc 



E 

/c=i 

iV 

E 

fe=i 



i#l 

8 Pfc 



Pfc + T^pfc - ^— ) ( 1 - T""; 



2 V 



fc=i 



8 Pfc 

A? 



fc=l 



8 Pfc 

JV 



1 N 

-^pudcj) 



fc=i 



1 W dp 2 N i N 1 N 

1 ~ a E ~T ~ 1 E PfcC ^ ~ 9 E ^fc^fc ~ o E Pk^l, 



8 Pfc - 
fc=i ^ K fc=i fc=i fe=i 

where we have made use of the normalization constraint and its differential form, 

N N 



J^Pfc = 1 and, ^ d Pk = 0, 



fc=i 



fc=i 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 
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that is, 



In (go (#)) d = 



dq 



(1 - qlY [(N 1) (1 - ql) + qlY 



de. 



(B2) 



The integral I (qo (9)) on the left-hand side of (|B2|) can be expressed in terms of an elliptic integral of the second kind 
E[-\-] (see [HI, for instance) as follows, 



In (qo (0)) = 



VJV (1 - gg) + 2gg - 1 
(N-l) 



(N — 2) q ~ ?o ^( Sin 1 (^o) l#zrf) 



N{ql-l)-2ql + l 



Ar(i-gg)+2gg-i 
JV-l 



(B3) 



Finding an analytical solution qo (9) satisfying the following equation, 

1 



In (go (0)) ~ 



Cn 



0. 



(B4) 



where Cn is a reaZ integration constant is quite challenging. Approximate solutions to Eq. ()B4|) would require 
numerical investigations which are left for elsewhere. However, in the two-dimensional case (N = 2), we get a closed 
form solution for qo (9). In such a case I2 (qo (9)) reads, 



X 2 (q (9)) = E (sin- 1 (q ) |0) d = arcsin (q (9)) , 



that is, 



q o (0) = sm(9 + C 2 ). 



(B5) 



(B6) 



As a conclusive remark, we stress that although finding an analytical solution qo (9) for Eq. (|B4|) may not be trivial, 
for sure go (9) in (IB6[) does not satisfy this equation when N ^ 2. This leaves open the possibility of finding actual 
parametric-dependent unitary probability paths also in the case of a single-item quantum search. 
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